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Abstract 

In this paper we report a theoretical model based on Green's functions, Floquet theory and averaging 
techniques up to second order that describes the dynamics of parametrically-driven oscillators with added 
thermal noise. Quantitative estimates for heating and quadrature thermal noise squeezing near and below 
the transition line of the first parametric instability zone of the oscillator are given. Furthermore, we give an 
intuitive explanation as to why heating and thermal squeezing occur. For small amplitudes of the parametric 
pump the Floquet multipliers are complex conjugate of each other with a constant magnitude. As the pump 
amplitude is increased past a threshold value in the stable zone near the first parametric instability, the 
two Floquet multipliers become real and have different magnitudes. This creates two different effective 
dissipation rates (one smaller and the other larger than the real dissipation rate) along the stable manifolds 
of the first-return Poincare map. We also show that the statistical average of the input power due to thermal 
noise is constant and independent of the pump amplitude and frequency. The combination of these effects 
cause most of heating and thermal squeezing. Very good agreement between analytical and numerical 
estimates of the thermal fluctuations is achieved. 

Keywords: parametric oscillator, parametric resonance, Floquet multipliers, Green's functions, averaging method, 
Langevin equation, and thermal squeezing. 
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I. INTRODUCTION 



Parametrically-driven systems and parametric resonance occur in many different physical sys- 
tems, ranging from the mechanical domain to the electronic, microwave, electromechanic, op- 
tomechanic, and quantum domains. In the mechanical domain we have Faraday waves [1], inverted 
pendulum stabilization, stability of boats, balloons, and parachutes [2]. A comprehensive review 
of applications in electronics and microwave cavities spanning from the early twentieth century up 
to 1960 can be found in Ref. [3]. A few relevant recent applications, in micro and nano systems, 
include quadrupole ion guides and ion traps [4], linear ion crystals in linear Paul traps designed 
as prototype systems for the implementation of quantum computing [5-7], magnetic resonance 
force microscopy [8], tapping-mode force microscopy [9], axially-loaded microelectromechanical 
systems (MEMS) [10], torsional MEMS [11]. In the quantum domain we could mention wideband 
superconducting parametric amplifiers [12] and squeezing in optomechanical cavities below the 
zero-point motion [13]. 

Parametric pumping has had many applications in the field of MEMS, which have been used 
primarily as accelerometers, for measuring small forces and as ultrasensitive mass detectors since 
the mid 80's [14]. An enhancement to the detection techniques in MEMS was developed by Rugar 
and Griitter [15] in the early 90's that uses mechanical parametric amplification (before transduc- 
tion) to improve the sensitivity of measurements. This amplification method works by driving the 
parametrically-driven resonator on the verge of parametric unstable zones. They were looking for 
means of reducing noise and increasing precision in a detector for gravitational waves, when they 
experimentally found classical thermomechanical quadrature squeezing, a phenomenon which is 
reminiscent of quantum squeezed states. The classical version is characterized by oscillating lev- 
els of the response of the parametric oscillator to noise at the frequency of the parametric pump, 
in such a way that the product between maximum and minimum output noise levels is constant. 
They observed that when the pump was turned on, the noise increased in one quadrature, while on 
the other it decreased. No theoretical model was proposed by them to explain the effect though. 
Subsequently, DiFilippo et al. [16] and Natarajan et al. [17] proposed theoretical explanations for 
this noise squeezing phenomenon, but their models did not treat noise directly in the equations of 
motion. 

Here, we study a parametrically-driven oscillator in the presence of noise with the objective of 
understanding what causes thermal squeezing and heating in the stable zone near the transition line 

2 



of the first parametric instability. The one-degree of freedom model studied here may be applied 
for instance to the fundamental mode of a doubly-clamped beam resonator that is axially loaded, 
in which case the one degree of freedom represents the amount of deflexion of the middle of the 
beam from the equilibrium position. The present model can also be applied to the linear response 
of ac driven nonlinear oscillators to noise (such as transversally-loaded beam resonators), see for 
example Ref. [18]. 

One of the objectives of the present investigation is to extend and improve on recently ob- 
tained analytical quantitative estimates of the amount of quadrature noise squeezing and heating 
in a parametrically-driven oscillator [19]. Here we use the Green's function approach, previously 
developed to solve the Langevin equation, aligned with averaging techniques up to second order, 
to obtain more precise analytical estimates of the thermal fluctuations in the parametrically-driven 
oscillator with added noise. We further show, using an approximate Floquet theory based on first 
and second-order averaging approximations, that thermal squeezing and heating are related to the 
onset of real- valued Floquet multipliers (FMs) with different magnitudes. It is shown that one FM 
grows while one gets closer in parameter space to the first transition line to instability while the 
other FM decreases. As a consequence, one gets two different effective dissipation rates, while at 
the same time the input power due to noise remains constant as the pump amplitude is increased. 
We show below that these effects account for most thermal squeezing and heating observed. Fur- 
thermore, first-order analytical estimates of heating and the amount of squeezing are also provided. 

The contents of this paper are organized as follows. In Sec. (II) we present our theoretical 
model, in Sec. (Ill) we present and discuss our numerical results, and in Sec. (IV) we draw our 
conclusions. 

H. THEORY 

The equation for the parametrically-driven oscillator (in dimensionless format) is given by the 
damped Matthieu's equation 

x + UqX = — 7X + F p cos(2u;t) x, (1) 

in which 7 and F p ~ 0(e), where e « 1. Since we want to apply the averaging method (AM) 
[20, 21] to situations in which we have detuning, it is convenient to rewrite Eq. (1) in a more 
appropriate form with the notation = luq — uo 2 , where we also have f2 ~ 0(e). With this 
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substitution we obtain x + u 2 x = —Qx — + F p cos(2ut) x. We then rewrite this equation in 
the form x = y,y = —u 2 x + f(x, y, t), where f(x, y, t) = —Vlx + F p cos(2ut) x — 7?/. We now 
set the above equation in slowly-varying form with the transformation to a slowly-varying frame 







[vl 





and obtain 




cos ut — sin ut \ I U 
-usmut — u cos ut I \ V 



cos cut — - sin tot \ l 
— smut — ^cosut I \ f(x,y,t) 

1 / sinutf(x,y,t) \ , , N / U 

" \cos utf{x, y,t) J \V 

The components of the jacobian matrix DF(U, V, t) of the above flow is given by 
DF n 
DF l2 
DF 21 



(2) 



(3) 



DF 



2u 
-1 
2u 
-1 
2u 
-1 



' ju(l - cos(2ut)) - sin(2wt) + y sin(4wt) 



22 



2u 



F 

Q(l - cos(2ut)) + -^(1 - 2 cos(2wt) + cos(4u;t)) + 70; sin(2u;t) 

+ cos(2u;t)) + -^[1 + 2cos(2o;t) + cos(4wt)] + 7cusin(2a;i) 
tt sin(2ut) + 70; (1 + cos(2wt)) - y sin(4wt) J 



(4) 



After application of the AM to first order (in which, basically, we filter out oscillating terms at 2u 
and Au in the above equation), we obtain 
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-1 
2u 
-1 
2u 



'yuu + ( + — - j r 



-VL + y ) u + •yuv 



(5) 



where the functions U(i) and V(i) are related to their slowly-varying averages u(t) and v(t), 
respectively, by the transformation 



(6) 



According to the averaging theorem [22, 23], the vector W obeys d t W(u, v, t) = g(u, v, t), where 
the vector g(u, v, t) corresponds to the explicitly-varying components of the right side of Eq. (3). 





1=1 


' u + Wi \ 






[v + W 2 / 



Namely, we have 



dW(u,v,t) 
dt 



gi(u,v,t) 
g 2 (u,v,t) 



1 / [(Q + F p )v + wyv] cos(2wt) + (Qu — ^yoov) sin(2o;t) — f[v cos(4u;t) + u sin(4o;t)] 
2uj I [(ft - F p )u - co-fv] cos(2wt) - (ttv + juu) sm(2ut) - \ [wcos(4wt) - vsin(4wt)] 

Upon integration we find 



W{u,v,t) 



4a; 2 



CJ7 sm(2wt) — Q cos(2cut) + cos(4o;t) 



(ft + F p ) sin(2c;i:) + 7 w cos(2wt) - ^ sin(4o;t) 



(ft — F p ) sin(2o;t) + cos(2wt) — ^ sin(4o;t) — ury sin(2wt) + ft cos^c^t) — ^ cos(4o;t) 



in which the integration constants are set to zero. The averaging theorem [21] states that these two 
sets of functions, namely (u(t),v(t)) and (U(t), V(t)), will be close to each other to order 0(e) 
during a time scale of 0(l/e) if they have initial conditions within an initial distance of 0(e). 
So by studying the simpler averaged system, one may obtain very accurate information about 
the corresponding more complex non-autonomous original system. Using the transformations 
u{t) = e-~ /t/2 u(t) and v(t) = e" 7 ' /2 'D(t) in Eqs. (5), we obtain 



u 



2uj V 2 



-1 

2lo 



2 



u. 



(7) 



Upon integration of Eqs. (7), one finds the solution 



u(t) 
v(t) 



-7*/2 



-7*/2 



1 / \ /3 — 5 , v 
mq cosh(Kt) H v smh(/CTj 



vq cosh(/tt) + 



K 

(3 + 6 



Uq sinh(Kt) 



(8) 



where k = \/ /3 2 — 5 2 , f3 = —F p /Auj, and 5 = Q/2oj. Hence, we find that the first parametric 
resonance, i.e. the boundary between the stable and unstable responses, is given by 



( 7W ) 2 = (Fp/2) 2 - ft 2 . 
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(9) 




)l 









This result is valid for oj rs uq even in the presence of added noise. In Fig. 1 we find very good 
agreement between the boundary obtained from numerical integration of Eq. (1) and the boundary 
given by the averaging technique. 

From Eqs. (2) and (8) we obtain the approximate fundamental matrix <3>(i) (also known as the 
time evolution operator) 

x{t) ) =*( f / Xo 

y(t) ) 

cos out — smut \ I cosh(/«t) — ^-sinh(/«t) 
-usmut —ucosut J \ ^-smh(Kt) — ^cosh(Kt) 

where $(0) = /, in which / is the identity matrix. From Floquet theory [20] we know that = 
P(t)e Bt , where P(t) is a periodic matrix with period T = n/co. We also know that -P(O) = /. 
The eigenvalues of e BT are known as the Floquet multipliers. We rewrite the fundamental matrix 
in the following form 

/ cos ut -smut\ I e" 7 * /2 cosh(«*) -^e~ 7i/2 sinhf/tt) \ 

$(t) « u \\ r ^ " K • (10) 

I -wsinwt coswt / I -^p e - 7i / 2 sinh(Kt) e" 7 */ 2 cosh(«t) / 

Hence, we notice that from the approximate solution of the fundamental matrix, via first-order 
averaging, we can find the approximate Floquet multipliers. They are given by 

A± = _ e -(i±«)7\ (U) 

Further improvements can be made by going to second order averaging. According to Ref. [23], 
the second-order corrections are given by the time average of 



DF(U,V,t)W(U,V) = 
- 1 I ~luF p 



- 7 V 



7 2 w 2 + n(n - f p ) + f 2 /8 



+ F p ) - F p 2 /8 



)l 









Hence, the second-order approximation to Eq. (3) becomes 



u 



-1 
-1 



/v 4u/ V 2 4 

7 2 , 



-0 + 



4w 2 



Q(Q + F p 



IT) 

F„ 



(12) 
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The solution is given by 

u(t) = e^ 2 
v(t) = e-^' 2 



, . sinh(££)\ ., . sinh(ft) 
cosh(et) + a u + (b - c) ^ 

. sinh(£t) / , . sinh(££)\ 
(b + c) + ( cosh(et) - a ^ J 



(13) 



with f = Va 2 + b 2 - c 2 , a = -g, b = f3 - 5/3/u, and c = 5 - 7 2 /(8u;) - 5 2 /(2w) - f3 2 /(4u). 
Thus, We can write the fundamental matrix in second-order approximation as 



$(£) f« e 



--ft/2 



cos cut -sin cut \ / cosh(^ t) + a sinh ^^ 



-^sinh(^) 



—cu sin cut cos cut 



sinh(£ t) cosh(£ t) — a- 



After some simple algebraic operations, we find the Floquet multipliers (eigenvalues of $(T)) to 
be given by 



A, 



,-(<y/2±0T 



(14) 



Hence, the transition line to instability in second-order averaging is given by 



7 = 2£. 



(15) 



A. Green's function method 



The equation for the Green's function of the parametrically-driven oscillator is given by 

d 2 



d_ 
Ft 



dt 2 + ^ + 7— -F p cos(2u;t) 



G(t,t') =5{t-t'). 



(16) 



Since we are interested in the stable zones of the parametric oscillator, for t < t' G(t, t') = and 
by integrating the above equation near t = t', we obtain the initial conditions when t = t' + + , 
G(t,t') = Oand §- t G(t,t') = 1.0. 



1. 1 st-order averaging 



Although Eq. (16) may be solved exactly by using Floquet theory [24], one obtains very com- 
plex solutions. Instead, we find fairly simple analytical approximations to the Green's functions 
and, subsequently, to the statistical averages of fluctuations using the averaging method. From 
Eq. (2) we obtain the approximate Green's function is G(t,t') = cos(ut)U(t) — sm(uit)V(t) = 
cos(ut)[u(t) +Wx{u(t),v(t),t)) - sxa(ut)[(v(t) + W 2 (u(t),v (£),£)], the functions u(t) andv(t) 
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are given by the solution of the system of coupled differential equations (8), where the time t is 
replaced by t — t' and the initial conditions set at t — t' are given by u(t') = — sm(cot')/uj and 
v(t') = — cos(ut')/u. For simplicity we set Wi(t') = W 2 (t') = 0. We find the approximate 
Green's function to be 



G(M') 



,-r(t-t')/2 



(Q ^ 
cos1i(k (t — t')) sm(ut') H sinh(/« (t — t')) cos( ut ; ) 
K 



sva(ujt) ( - sinh(/t (t — t')) sin(wt / ) + cosh(/t (t — t')) cos(ut'] 



- 7 (t-f)/2 



5 



cosh(K (t - f)) sm(u(t - t')) + - sinh(« (t - t')) cos(w(t - i 7 )) 
a; ^ K 

— — sinh(/« (t - t')) cos( u(t + £')) 

K 



(17) 



for t > t' and G(t, t') = for t < t' . In the stable zone of the parametrically-driven oscillator, 
when \(3\ > \5\ , we can rewrite the Green's function replacing the initial conditions and using 
simplifying trigonometrical identities. The change of variables t' — t — r leads to 



e" 7r/2 f 5 
G(t, t — t) ps < cosh(/t r) sin(a;r) H — sinh(/« r) cos(wr) 

W t K 

— — sinh(K r) cos(u;(2t — r)) 



(18) 



2. 2nd-order averaging 



The Green's function obtained by second-order averaging is given by 



- 7 (t-t')/2 

G(t,0 « x 

CO 



, m smh(£(t-t'))\ • , a /, xSinh(^(t-t')) , a' 
cosh(f (t - t')) + a ^- sm(ut') + (b - c) ^- - cos(W) 



cos(cut) 

'Shi(wt) , ( , ^ ,, 

for t > t' and G(t, t') = for t < f. This can be rewritten in a shorter format as 



,sinh(f(t-f)) . , A / , m sinh(C(i-0)'i 
(6 + c) ^- sm(wt') + cosh(£(t - £')) - a ^- - cos(a)i') 



e - 7 r/2 r sinh(^r) 
t — r) ps < cosh(£r) sin(wr) + c : cosier ) 



sinh(^r) 



[asin(w(2t - r)) + bcos(cu(2t - r)] 



(19) 



for t > and G(t, t — r) = for r < 0. We notice that this second-order approximate Green's 

function can be put in the same format as the one of the first-order approximation given by Eq. (18). 

„, . e~ 7r / 2 f . N . . sinh(£r) . . 

Git, t — t) ps < cosnffr) smiuT) + c cos(wr) 

u { £ 

-/?' S -^^l cos(u[2(t - t ) - r]) } , (20) 



where (3' = (3^(1 - S/tu) 2 + 7 2 /(4w 2 ), a = sin(2wt ), and 6 = /?' cos(2wt )- In Eq. (18) we 
have to replace k by ^, 5 by c, /3 by and t by t — to to obtain the second-order Green's function 
given by Eq. (20). Note that the Green's functions given here represent the first two steps of a 
Green's function renormalization procedure based on the averaging method. 

B. Thermal fluctuations 

We will now investigate the effect of noise on the parametric oscillator [19]. We start by adding 
noise to Eq. (1) and obtain 

x = — (jJqX — 7± + F p cos(2o;t) x + R(t), (21) 

where R(t) is a random function that satisfies the statistical averages (R(t)) = and (R(t)R(t')) = 
2T^S(t — t'), according to the fluctuation-dissipation theorem [25]. The temperature of the heat 
bath in which the oscillator (or resonator) is embedded i T. Once we integrate these equations 
of motion we can show how classical mechanical noise squeezing and heating occur. We shall 
now review the analytical method developed in Refs. [19, 26] to study the parametrically-driven 
oscillator with noise as given by Eq. (21). 

Using the Green's function we obtain the solution x(t) of Eq. (21) in the presence of noise R(i) 

/oo 
dt' G(t,t')R(t'), (22) 
■oo 

f°° d 

x(t) = v(t) = v h (t) + J dt' -G{t, t')R(t'), (23) 

where Xh(t) is the homogeneous solution, which in the stable zone decays exponentially with 
time; since we assume the pump has been turned on for a long time, Xh(t) = 0. By statistically 
averaging the fluctuations of x(t) we obtain 

r poo poo 

( x 2 (t))= dt' dt"G(t,t')G(t,t")(R(t')R(t")) =2T 7 / drG(t,t-T) 2 , (24) 

JJ-oo JO 



-oo 

ft 



(v 2 {t)) = 2T 7 / dt' 



oo 



(25) 



where r = t — t' . 

We notice that by varying the pump amplitude F p and the detuning Q, we can create a contin- 
uous family of classical thermo-mechanical squeezed states, generalizing the experimental results 
of Rugar and Griitter [15]. An estimate of the time average of the statistically averaged thermal 
fluctuations, when \/3\ > \5\,is given by 



<**(*)) 



2T 7 



./o 



-7r 



cosIi(kt) sin(wr) H — sinh(Kr) cos(a;r 



2 ^ 
f 2k 



^ [h + I 2 + h + h\ 



sinh 2 (kt) ^ dr 



(26) 



where the integrals are given by 



h 
h 



2k 2 



e 7T sinh 2 (kt) fir 



/? 2 







7(7 2 — 4k 2 ) 



e 7r cosh (kt) sin (u)t)cIt 


1 f 1 7 

+ 



7 



1 



2 [2 7 2(7 2 + 4w 2 ) 2( 7 2 -4k 2 ) 4 

^2 /-co 



Re 



1 



1 



K' 



/■oo 

/ e~ JT sinh 2 (Kr) cos 2 (uT)dr 
Jo 



1 



1 



+ 



7 — 2k — 2%uj 7 + 2k — 2iio 



7 



k 2 [ 7(7 2 — 4k 2 ) 

S f°° 1 
— / e~ 7T sinh(2Kr) sin(2u;r)c?T = -Im 

2k ./n 



7 — 2k — 2zo; 7 + 2k — 2io; 

1 

Im 

4 



4( 7 2 + 4w 2 



k(7 — 2k — 2zw) k(7 + 2k — 2iu) 



An estimate of the statistically averaged thermal fluctuations, when \/3\ > \5\, is given by 



(x 2 (t)) « (x 2 (t)) + A 2w cos(2o;i:) + B 2u] sm(2ut) + A 4uJ cos{4ojt) + B 4uj sin(4o;t), (27) 



where 



A 2uJ = -^(K 1 + K 2 ) 
B 2u) = -^(K 3 + K,) 
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with 



K 3 
K 4 



Im 



5 



1 

k(7 — 2k, — 2iu) k(7 + 2k — 2iuj) 
7 1 



K 2 



Re 



1 



4k 



Sk 
8k 2 
4k 2 



( 7 2 - 4k 2 ) 



Re 



1 



7 — 2k — 2iu 7 + 2k — 2iu 
1 



-Im 



7 — 2k — 2iw 7 + 2k — 2icu 
1 2 



7 — 2k — 2%uj 7 + 2k — 2%u 7 — 2icu 

1 1 2 

+ 



_(7 + 2k) 2 + 4w 2 ( 7 -2k) 2 + 4u; 2 7 2 + 4w 2 _ 
The remaining coefficients of Eq. (27) are given by 

/3 2 T7„ r 1 1 

A 4uJ = ^^Re 



B 



4tU 2 K 2 

/3 2 T 7 



4w 



4wk 2 



Im 



7 — 2k — 2iw 7 + 2k — 2iu 7 — 2ico 
1 1 2 



7 — 2k — 2iw 7 + 2k — 2iu; 7 — 2iu 

Notice that when one gets close to the zone of instability, we obtain a far simpler expression for 
the average fluctuations. It is given approximately by 



<* 2 (t)> 



2T 



2 ' 7 + * 2 



/3 2 + 



4/3T 



7 2 — 4k 2 cj 2 (7 2 — 4k 2 ) 
It is easy to verify that the minimum of the position fluctuation is given by 

T v/<5 2 + 7 2 /4-|/3| 



7 

5 cos(2wt) + — sin(2o;t) 



a. 



min 



= (a; 2 (t)> 



min 



and the maximum is given by 



2 

°max 



(x 2 (t)) 



max 



such that 



2w 2 ^2 + 7 2/ 4+ ^1 



r y/^ 2 + 7 2 /4+|/3| 
2w 2 ^2 + ^/4 _ 



a min a max ~ 



(28) 



(29) 



This is a verification that the classical phenomenon of thermal squeezing indeed occurs near the 
transition line to instability. These approximate expressions for squeezing are not valid with de- 
tuning (5 ^ 0) when ,5 = 0, although they correctly predict no squeezing in such situation and the 
equipartition theorem is correct to 0(e) if VI = 0(e). 
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From Eqs. (20) and (24) we obtain the dc contribution to fluctuations in the second-order 
approximation 



(x 2 (t)} « ^ f dr <~ 7 



i )2 

U j 



cosh(£r) sin(wr) + c 



sinh(^r) 



cos(cjt) 



, ^ 2 sinh 2 (Cr) 



2£ 2 



2T 7 



/l + / 2 + / 3 + ^4 



(30) 



where the integrals are given by 

/3 12 



h 
h 



7(7 2 -40' 
1 f 1 7 



+ 



7 



2 [2 7 2(7 2 + 4cu 2 ) 2( 7 2 -4£ 2 ) 4 



Re 



+ 



1 



e 



i 



e b(7 2 -4e) s 



+ 



7 - 2£ - 2iu 7 + 2£ - 2icu 



7 



/ 4 = -Im 



7 - 2£ - liuj 7 + 2£ - 2iw 
1 



4(7 2 + Au 2 ) 



,e(7 - 2£ - 2zu;) £( 7 + 2£ - 2iw) 



All the other first-order approximation results given at Eqs. (27-29) are the same in second-order 
approximation, except for the parameter replacements given at the end of subsubsection (II A 2). 
It is noteworthy to mention that the squeezing condition given at first-order approximation in 
Eq. (29) is still valid in the second-order approximation. This implies that a renormalization 
procedure based on the averaging method is possible for the Green's function of the parametric 
oscillator with parameters set near the onset of the first instability zone. 



C. Energy balance 

It is important to verify that in the stable zone of the parametric oscillator, on average, the power 
input from the external noise and from the internal pump is balanced by the dissipated power. 
The instantaneous input power due to the additive noise is given by 

POO Q 

Pnoiseit) = x(t)R(t) = v h (t)R(t) + J dt' —G (t, t') R(t') R(t) . (31) 
Hence, we obtain the statistically averaged noise input power 



(Pnoiseit)) = J dt' ^G(t,t')(R(t')R(t)) = 7 T 



= lT. (32) 

t'=t- 
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The statistically-averaged dissipated power is given by 



(P dlss (t)) = -l(m 2 ) = -l I dU I dt" ^G(t,t')^G(t,t")(R(t')R(t")) 



oo J — oo 
2 



= -2 7 2 T f dt' 

J —oo 

The statistically-averaged pump power is given by 



(33) 



(Ppum P (t)) = F p cos(2ut)(x(t)x(t)) 

= F p cos(2ujt) j dt 1 j dt" G(t,t')^-G(t,t")(R(t')R(t")) 

J -oo J -oo 

= 2 7 TF p cos(2u;i) / dt' G(t,t')^-G(t,t') 

J-oo 

= 1 TF p cos{2ut)^- f dt' G(t,t') 2 

J —oo 

= ^cos(2^)|(x(t) 2 >. (34) 

From Eq. (27) we obtain approximately the time-averaged statistically-averaged pump power. It 
is given by 

(Ppump(t)) ~ — — E>2u)- (35) 

Near the threshold to instability in first-order approximation it becomes 

F 2 T 7 

(Ppum P (t)) « 4w 2( 7 2_ 4K 2)' 

while on second-order approximation it is 



In the stable zone of the parametric oscillator, when the stationary point is reached, one gets 
the energy balance on average, that is 



{Pnoiseit)) + (P diss (t)) + (P„(t)> = 0. (36) 

One can use the above expression to obtain the time average of the velocity fluctuations 

wm = (t+ (wi/t) = (t+ ^bJj . 

By comparison of the above equation with Eq. (30) one notices that the equipartion of energy 
breaks down once pumping is on. 
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HI. RESULTS AND DISCUSSION 



In Figs. (2, 3, 4) we show the dependance of the magnitude of the Floquet multipliers (FMs) 
on the pump amplitude F p . When the FMs given by Eq. (14) become real, they branch off in two 
magnitudes. When this occurs one of the FMs (A_) becomes larger as F p is increased, eventually 
becoming larger than one, when the system given by Eq. (1) becomes unstable, while the other FM 
(A+) becomes smaller. From Eq. (14) this implies into two effective dissipation rates, along the sta- 
ble manifolds of the Poincare first-return map (with period n/co, i.e. $(7r/u;)). This phenomenon 
causes both heating and squeezing of the parametric oscillator with added noise, since the average 
power input due to thermal noise remains constant as is shown in Eq. (32) and one effective dis- 
sipation rate is decreased. It also causes quadrature thermal squeezing since in one direction, the 
A_-stable manifold of $(n/u)) in the phase space of u and v, there is less effective dissipation and, 
consequently, more fluctuations, while in another direction, the A + -stable manifold of $(ir/u) in 
the phase space of u and v, there is more effective dissipation and , consequently, less fluctuations. 
This inbalance in the effective dissipation rates, we claim, is the main cause of thermal squeezing 
in the parametric oscillator with added noise. The dependance on pump amplitude of the effective 
dissipations can be seen on Fig. 5 (second-order averaging result) and on Fig. 6 (Floquet theory 
numerical result). 

In Fig. 7 we show several squared Green's functions with initial conditions spread out evenly in 
time during one period of the pump (ir/cu). They are vertically spaced only for clarity, since all of 
their assymptotes are zero. The Green's functions are shown squared because that is the way they 
contribute to the thermal fluctuations in Eq. (24). One notices that the second-order approximation 
Green's functions yield a much better approximation to the numerical Green's functions than the 
first-order approximation. This is specially evidenced the closer one gets to first transition line to 
instability, a consequence of the fact that the second-order analytical expression for the transition 
line given in Eq. (15) is a better approximation than the first-order expression in Eq. (9). 

In Fig. 8 we show a comparison between time series of the fluctuation (x 2 (t)) given by Eq. 
(24) in which the Green's functions are given either by the second-order approximation expres- 
sion from Eq. (19) or by the numerical Floquet theory Green's functions. We use several different 
Green's functions with negative detuning (a; = 0.95), in resonance (u = 1.0) and positive detun- 
ing (oj = 1.05) all with the same pump amplitude F p = 0.185. We observe that the numerical 
and approximate Green's functions are very similar and that the squeezing amplitude and heat- 
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ing (proportional to the time average of (x 2 (t))) are very dependant on detuning from resonance. 
From Fig. 1 one sees that the off resonance results are slightly below the threshold of the strong 
heating and squeezing zone. In Fig. 9 we show again several time series of the fluctuation (x 2 (t)), 
but this time the pump amplitudes are chosen such that the parameters are inside the heating zone 
(as given in Fig. 1) and close to the transition line to instability. One sees then considerably higher 
dynamical temperatures and squeezing amplitudes for the detuned fluctuations than in the previous 
figure. 

In Fig. 10 we show a logarithmic plot of the dc component of the mean-square displacement 
over the heat bath temperature. Most of the heating occurs inside the heating zone in which the 
Floquet multipliers are real, each one with a different amplitude, one that increases and the other 
that decreases as the pump amplitude is increased. 

In Fig. 1 1 we show level sets in decibels of the dc component of the mean- square displacement 
over the heat bath temperature in second-order approximation. One sees that most of the heating 
occurs inside the heating zone in which the Floquet multipliers are real. In Fig. 12 we show level 
sets of the squeezing amplitude of the mean-square displacement over the heat bath temperature in 
second-order approximation. This attests that most of the squeezing also occurs inside the heating 
zone as claimed before. 

IV. CONCLUSIONS 

Here, we studied a parametrically-driven oscillator with added noise with the objective of un- 
derstanding what causes heating and thermal squeezing in the stable zone near the transition line 
of the first parametric instability. We improved on our previous work [19] and obtained a more 
accurate expression for the Green's functions and, consequently, obtained more precise estimates 
of the amount of heating and quadrature noise squeezing in the parametrically-driven oscillator. 
Furthermore, we used an approximate Floquet theory, based on first and second-order averaging 
approximations, to explain why heating and thermal squeezing occur in the parametrically-driven 
oscillator with added noise investigated here. These phenomena are related to the onset of real- 
valued Floquet multipliers (FMs) with different magnitudes. It was shown that as one FM grows 
while one gets closer (in parameter space) to the first transition line to instability the other FM 
decreases. As a consequence, one gets two different effective dissipation rates, while at the same 
time the input power due to noise remains constant as the pump amplitude is increased. We showed 
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that these effects account for most thermal squeezing and heating observed. We showed as well 
that the second-order Green's function of the parametric oscillator has the same form as the first- 
order Green's function, which implies that with a simple parameter change given in the text all 
first-order results for the amount of heating and thermal squeezing also apply, with increased ac- 
curacy, to second-order approximation. This indicates that a Green's function renormalization is 
very well feasible if one goes to higher orders of approximation in the averaging method. 

The one-degree of freedom model studied here may be applied for instance to the dynamics of 
the fundamental mode of an axially loaded doubly-clamped beam resonator. The present model 
could also be applied to the linear response of ac driven nonlinear oscillators, such as transversally- 
loaded beam resonators, to noise. We note further that this model can be applied as well to the 
investigation of the dynamics of ions in quadrupole RF ion guides or traps [4] in the presence of 
thermal noise. Most importantly, the notion that we can have two different effective dissipation 
rates, one for each stable Floquet multiplier, could be used to create control schemes similar to 
those of the OGY method [27] to reduce the effects of noise even further in parametric amplifiers 
and in RF ion traps. 
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FIG. 1. (color online) Comparison between numerical and averaging method predictions for the boundary 
of the first instability zone of the damped parametrically-driven oscillator of Eq. (1). In the region above the 
thick solid black line lies the unstable zone obtained by numerical computation, the red dashed line is the 
second-order averaging prediction [Eq. (15)] for the transition line. The numerical results are obtained by 
numerically calculating the corresponding Floquet multipliers. The transition line to parametric instability 
is defined when at least one of them has modulus equal to 1. The heating zone occurs when the Floquet 
exponents given in Eq. (14) become real. The heating boundary (blue long dashed-dotted line) is equivalent 
to the transition line to parametric instability when 7 = 0. 
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FIG. 2. (color online) The magnitude of the Floquet multipliers at ui = 0.9,7 = 0.1. The transition to 
instability of the dynamics of the parametric oscillator occurs when one of the Floquet multipliers becomes 
bigger than 1. The second-order averaging prediction [Eq. (14)] of the FMs reproduces very accurately the 
numerical results. The branching off of the FMs occurs when they become real. 
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Magnitude of Floquet Multipliers oo=1.0 




FIG. 3. (color online) The magnitude of the Floquet multipliers at cj = 1.0, 7 = 0.1. From the onset of the 
parametric drive the FMs are real. The second-order averaging prediction [Eq. (14)] of the FMs reproduces 
very accurately the numerical results. 
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FIG. 4. (color online) The magnitude of the Floquet multipliers at oj = 1.1, 7 = 0.1. The second-order aver- 
aging prediction [Eq. (14)] of the FMs reproduces very accurately the numerical results. Off resonance one 
notices that the pump amplitude has to be higher before the branching off of the Poincare map eigenvalues 
(FMs) in agreement with the results of Fig. (1). 
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Effective dissipation rates 




FIG. 5. (color online) Effective dissipation rates based on second-order averaging approximation obtained 
from Eq. (14). Heating and squeezing mostly occur when the floquet multipliers become real and branch 
off in two different values. Consequently, there is an effective dissipation rate associated with each floquet 
multiplier. An effective dissipation rate 7 e jj < 7 cause heating. The different values of j e ff, one smaller 
than 7 and the other larger than 7, result in squeezing. 7 = 0.1 
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Effective dissipation rates 
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FIG. 6. (color online) Numerical effective dissipation rates. Same parameters as in Fig. (5). 
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FIG. 7. (color online) In the frames above we show several squared Green's functions with equally-spaced 
in time initial conditions in one given period of the parametric driving. They are vertically spaced for clarity, 
since all asymptotes are zero. In each frame we have a comparison between numerical results given by the 
numerical integration of Eq. (16) and the analytical approximate results given by Eqs. (17) or (19). We have 
a) en = 0.9, b) uj = 1.0, and c) w = 1.1. The initial values of the Green's functions are G(t, t') = and 
§- t G(t, if) = 1.0 when t = t? + 0+ The pump amplitudes are indicated in the figure frames. 
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mean square deviation, 7=0.1, F p =0.185, T=1.0 
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FIG. 8. Mean-square deviation {x 2 (t)} time evolution. Comparison between Floquet theory numerical 
results and second-order averaging results given by Eq. (24) with the Green's function given by Eq. (19). 
The dot-long dashed line indicates the temperature of the external heat bath. All pump amplitudes are the 
same F p = 0.185. This indicates that the close one is to resonance the larger the squeezing amplitude and 
average dynamical temperature. 
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FIG. 9. Mean-square deviation {x 2 (t)} time evolution. Comparison between Floquet theory numerical 
results and second-order averaging results given by Eq. (24) with the Green's function given by Eq. (19). 
The dot-long dashed line indicates the temperature of the external heat bath. The values of the pump 
amplitude were chosen to be close to the transition line given by Eq. (15) 
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FIG. 10. (color online) Log plot of the dc component of the mean square displacement in second-order 
approximation, (x 2 (t)}, as given by Eqs. (30). In the parametric oscillator with thermal noise the dynamic 
temperature of the oscillator grows monotonically until it diverges at the transition line between stable and 
unstable zones. The simplified approximating curves are given by Eq. (28). 
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FIG. 11. (color online) Countour log plot of dynamic temperatures (the dc component of the statistical 
average of the square displacement in second-order approximation, {x 2 (t)} , as given by Eqs. (30)) over 
heat bath temperature. Most of the heating occurs inside the heating zone as predicted by our hypothe- 
sis. Each level curve is a dynamic isothermal. The levels are given in decibels given by the expression 
10 log((x 2 (t)} /T), hence each increment of approximately 3 dB increases the dynamic temperature by a 
factor of 2. 
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FIG. 12. (color online) Countour log plot of the squeezing amplitude (the dc component of the statistical 
average of the square displacement in second-order approximation, \J A\ w + B\ w , as given by Eqs. (27)) 
(with the appropriate second-order corrections) over heat bath temperature. Most of the squeezing occurs 
inside the heating zone as predicted by our hypothesis. Each level curve has constant squeezing amplitude. 
The levels are given in decibels given by the expression 10 log(yM| w + B^/T). 
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